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Spectral index properties of milli Jansky radio sources 
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ABSTRACT 

At the faintest radio flux densities (Si. 4 < lOmJy), conflicting results have arisen regard- 
ing whether there is a flattening of the average spectral index between a low radio frequency 
(325 or 610MHz), and e.g. 1.4GHz. We present a new catalogue of 843 MHz radio sources 
in the ELAIS-S1 field that contains the sources, their ATLAS counterparts, and the spectral 
index distributions of the sources as a function of flux density. We do not find any statistically 
significant evidence for a trend towards flatter spectral indices with decreasing flux density. 
We then investigate the spectral index distribution with redshift for those sources with reli- 
able redshifts and explore the infrared properties. An initial sample of faint Compact Steep 
Spectrum sources in ATLAS is also presented, with a brief overview of their properties. 
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1 INTRODUCTION 

The desire to understand the properties of the faintest radio source 
populations has led to numerous surveys pushing to ever fainter 
levels. There are now a large number of deep and wide area 
radio surveys availabl e, such as the Australia Telescope Larg e 
Area Survey (ATLAS ; lNorris et al.ll2006l : lMiddelberg et alj|2008?) 
the Cosmic Evolut ion Survey (COSM o"sTlScoville et afl 120071; 
Smolcic et ai]|2008l) . the ATESP Survey dPrandoni et alj|2000 jlbl 



2006), the Phoenix Deep Survey ( Hopkins et aL 2003), and many 



other s dSevmour et ail 2008; Ow en & Morrisonll200i : lOwen et al.l 
l2009t Irbar et alJl2(XH |201Ch . The bright radio source population 
(Si. 4 > lOmJy) is well studi ed, and is predominantly composed of 
active galactic nuclei (AG N) (Condon 1984; Gruppioni et al. 1999 



Maglio cchetti et al.l |2000| : Georgakakis et alj 1 1 9991 : lAfonso et al.l 
2006). At fainter flux densities, star-forming galaxies (SFGs) begin 



to domi nate the radio source population, particularly into the /xJy 
regime dwindhorst et ai][l9 85; Hopki ns et ail 1 19981: lAfonso et all 
l2005l ; ls"evmour et alj2008h . Understanding these faint radio source 
populations is essential for understanding galaxy evolution, the role 
of star formation, AGN and the relationship between the two. The 
use of multiwavelength data to complement radio surveys is vital in 
order to distinguish between AGN and star formation processes as 
the origin of the radio emission down to the faintest flux densities, 
although there is likely also a composite population at these faint 



levels dHill et al"I l999. 2001). The relative proportions of these two 
populations will have an effect on the average radio spectral index 
(a]j as a function of flux density, and at the faintest flux densities, 
core-dominated AGN may be more prevalent, and may flatten the 
average spectral index to a > —0.7. 

Conflicting evidence has arisen over the nature and proper- 
ties of sources at frequencies below 1.4 GHz, particularly of their 
spectral index properties, and whether there is a flattening of the 
average spectral indices for faint (5*1.4 < lOmJy) radio sources. 
IPrandoni et ail d200d [201 l|) found sources with fluxes less than 
a few millijansky had an average spectral index which was flat- 
ter than that of the brig hter radio sources (a ~ —0.7). Similarly, 
lOwen & Morrison! d2008TJ found a flattening of the average spectral 
index between 325 MHz and 1.4 GHz for 1.4 GHz selected radio 
sources with Si. 4 < lOmJy, and angular sizes > 3". They found 
though, that the spectral indices steepen again at the faintest flux 
density end of the 20 cm survey (~ 0.5 mjy). In con trast, in the 
deepest radio field, the Lockman Hole dlbar et al . 2009), no flatten- 
ing of the spectral indices between 610 MHz and 1.4 GHz was seen 
for flux densities Si. 4 > 100/Hy. 

A flattening of the average spectral indices implies that there 
is a flat or inverted spectrum population of sources at these milli- 
jansky flux densities, at fluxes fainter than where the steep spec- 
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1 Assuming Si, oc u a , where S is the measured flux density and v is the 
observer's frame frequency. 
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tram star-forming population emerges (~ 0.5 mjy). Investigating 
the millijansky and microjansky radio source populations at dif- 
ferent frequencies allows us to investigate this suspected flattening 
of the average spectral indices, and the population of sources caus- 
ing the flattening, particularly if we can fill in the frequency regime 
from very low frequencies (~ 100 MHz) up to 1.4 GHz. It is nec- 
essary to understand the spectral index properties of these faint ra- 
dio sources for many reasons, such as the z-a relation, used to find 
the most distant radio galaxies dde Breuck et alj20 04; Klam eTet al.l 
l200dllshwara-Chandra et alj2010h , identifying young radio AGN, 
such as Gigahertz Peake d Spectrum (GPS) and Compact Steep 
Spectrum (CSS) sou rces dO'Dedl 19981 : iPolatidis & Conwavll2003l : 
Mor ganti etaiT2 009) or determining the emission mechanism in ra- 
dio galaxies, whether it is from star formation, or the central AGN 
jPrandoni et alj|2006h . 

Here we present radio observations at 843 MHz in the ELAIS- 
Sl field, for which we have complementary 1.4 and 2.3 GHz data 
from ATLAS that we can use to investigate the possible spectral 
index flattening at a frequency nearer to the ubiquitous 1 .4 GHz ra- 
dio surveys. Section|2]briefly reviews the existing multiwavelength 
data covering the ELAIS-S1 region and ATLAS, and in Section [3] 
we describe our observations, the data reduction, and the cross- 
matching process to the ATLAS 1.4 and 2.3 GHz catalogue. We 
present the catalogue in Section [4] Our results and analysis are 
explored in Section [5] and a new faint sample of candidate CSS 
sources is presented in Section [6] discussing the initial selection 
and properties of these sources. Our results and plans for future 
work are discussed in Section|7]and our conclusions are presented 
in Section[8] Throughout this analysis, unless otherwise noted, we 
use the cosmological pa rameters, Qm = 0- 27, SI a = 0.73 and 
H = 71kms" 1 Mpc" 1 jKomatsu et ai]|201lh . 



2 THE ELAIS-S1 REGION 

The European Large Area ISO Survey - South 1 Region (ELAIS- 
Sl) has been a target for many multiwavelength surveys and tar- 
geted observations over the last decade. The ELAIS-S 1 field covers 
2° x 2°, centred on RA = 00 h 34 m 44 s .4 and Dec=-43°28'12".0 
(J2000.0). Observations covering the ELAIS-S 1 field include ra- 
dio imaging, optical imaging and spectroscopic redshifts, infra-red 
observations (near, mid and far-infrared), UV, and X-ray observa- 
tions. Part of the attraction for the multiwavelength surveys in this 
region is that this field ha s the lowest Galactic 100/xm cirrus e mis- 
sion in the southern sky l iSchlegel. Finkbeiner & Davisll 19981) . in- 
cluding the absolute minimum. 

2.1 Multiwavelength data in ELAIS 

2.1.1 ISO 

The field was first observed as part of a deep, wide-angle sur- 
vey with the Infr ared Space Observatory (ISO) at 15 and 90 /im 
dOliver et al.|[200(il) . Follow-up observations were done at 6.7 /im, 
and the catalogue cov ering the SI Region at these three wave- 
lengths is discussed in IRowan-Robinson et al.l fe004l) . The com- 
plete 15 pirn catalogue contains all the sources fro m the ELAIS 
region s; Nl, N2, N3, SI, S2, and was finalized by IVaccari et al.l 
(2005). The three bands of this survey have rms noise levels of 1.0, 
0.7 and 70 mjy respectively for the 6.7, 15 and 90 /im bands. Photo- 
metric uncertainties for all bands were ~ 10%, and the astrometric 
accuracy is ~ 0.5". 



2.1.2 Spitzer/SWIRE 

The ELAIS-S 1 field was observed by the Spitzer Space Telescope 
as part of the largest Legacy Progra m, the Spitzer Wide-area Infra - 
Red Extragalactic survey (SWIRE; lLonsdale et ail 120031 , l2004t) . 
SWIRE aimed to trace the evolution of dusty SFGs, AGN, and 
evolved stellar populations out to a redshift of z ~3, by imaging 
large areas of sky in seven different infrared bands. In combination 
with optical imaging, this would allow SED modeling and explo- 
ration of galaxy evolution with environment. The seven imaging 
bands of SWIRE used two different instr uments aboard Sp itzer, 
the Multiband Imaging Photometer (MIPS; Riek e et al|2004l) . and 
the Infrared Array Camera (IRAC; iFazio et alj 120041) , covering 
~ 7 deg 2 over the ELAIS-S 1 region. The 24 /im MIPS imaging 
reaches a 5a sensitivity of 350 /ily, and the IRAC bands of 3.6, 4.5, 
5.8 and 8.0 fim have 5a sensitivities of 3.7, 5.3, 48 and 37.75 ^iJy 
respectively. 



2.1.3 ESIS 

The ESO-Spitzer Imaging extragalactic Survey (ESIS) is the opti- 
cal follow-up to SWIRE, consisting of optical imaging in the B, V, 
R bands with the Wide Field Imag er (WFI;lBerta et al J2006!) , and I 
and z band imaging with VIMOS ( iBerta et al.l2 008). Currently, the 
VIMOS observations have been completed, but the WFI observa- 
ti ons are not yet fully processed. The BVR observations described 
in Bertaet all d2006T) cover 1.5 deg 2 , in the central region of the 
ELAIS-S 1 field. The catalogue of BVR sources is 95% complete 
to 25 m in B and V, and 24.5" 1 in R. 132712 sources are included 
in this catalogue, with an rms uncertainty of ~ 0.15" for the co- 
ordinates of the sources. The VIMOS data covers ~ 4 deg 2 in the 
I band and ~ 1 deg 2 in the z band, resulting in a completeness of 
90% at 23. l m in the I band, and 22.5 m in the z band. Over 300,000 
sources were catalogued in the I band, and over 50,000 in the z 
band, with an rms of ~ 0.2" in both bands. 



2.1.4 Previous Radio Observations 

Prior to the ATLAS observations (detailed below in §2.21 , 
iGruppioni etahl dl999L G99) observed the ELAIS-S 1 field using 
the Australia Telescope Compact Array (ATCA) in 1997, covering 
~ 7 deg 2 . The observations consisted of a mosaic of 49 different 
pointings, resulting in an image rms of ~ 80/iJy at 1.4 GHz, and 
a catalogue of 581 radio sources. The lowest flux density sources 
were catalogued down to a 5a level of 0.2 mjy in the centre of the 
field, and 0.4 mjy in the remaining area. 



2.1.5 X-ray Observations 

BeppoSAX, an X-ray satellite, first o bserved 40% of the ELAIS-S 1 
field in 1999 jAlexander et alj200ll) . These observations reached a 
sensitivity of ~ 10~ 13 ergcm _2 s _1 in the 2— lOkeV range, cov- 
ering ~ 1.7 deg 2 . The central ~ 0.6 deg 2 were then surveye d 
by XMM-Newton in four deep pointings JPuccetti et alj 12006). 
The XMM observations were taken in both soft and hard X-rays, 
and each pointing had a net exposure time of ~ 60 ks. A total 
of 478 sources were detected, with 395 in the soft X-ray band 
(0.5-2keV) and 205 in hard X-ray (2-10keV). The flux limits are 
~ 5.5 x 10" 16 ergcm" 2 s _1 and ~ 2 x 10" 15 ergcm" 2 s _1 re- 
spectively for the soft and hard bands. 
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Figure 1. The greyscale MOST 843 MHz image, with overlays indicating the approximate borders of the ATCA 1 .4 GHz and 2.3 GHz observations. The larger, 
red boundary represents the 1.4 GHz mosaic. 



2.1.6 UV Observations 

The u ltraviolet Galaxy Evolution Explorer (GALEX; lMartin et al.l 
120051) observed the ELAIS field in two of its observing 
modes; the Deep Imag ing and Wide Spectroscopic Survey modes 
(Burg areila et al.l ; 2005). The bands observed were centred on 153 
and 23 1 nm, and the current data release covers 70% of 

the sky in the imaging mode, and has 61439 spectroscopic sources. 
The ELAIS field was observed for ~ 9 hours in the spectroscopic 
survey mode, and ~ 3 hours in the imaging mode. 

In the analysis below we focus on the infrared data as the pri- 
mary complement to our radio data, although in future work we 
will also take advantage of the X-ray and UV measurements. 



2.2 The Australia Telescope Large Area Survey (ATLAS) 

ATLAS is the widest, d eep radio survey to date dNorris et al.l2006j : 
Middelber g et alj I2008L N06, M08), covering ~ 7deg 2 over two 
fields, ELAIS-S1, and the Chandra Deep Field South (CDFS). The 
rms of the 1.4 GHz imaging data is currently 30 /xJ y, and the aim 
is to ac hieve an rms of 10 /xJy across both fields (Banfie ld et al.l 
lin prepj) . ATLAS observations have been completed at 1.4 and 2.3 
GHz with ATCA from 2006 to 2010, with one data rele ase in 2008, 
and subsequent data releases to begin in late 2011 (Hal es et"al] 



http://galex.stsci.edu/GR6/ 



in prep.; Banfielc Tet aljin prepj) . ATLAS has many scientific goals, 
primarily focussed on investigating the evolution of galaxies and 
AGN. Specific goals include distinguishing between AGN and 
SFGs and determining the contribution of each to a given galaxy's 
luminosity, searching for high- z radio galaxies t o trace the forma- 
tion of clusters at high redshift (Mao et aL 2010), and fin ding new 
types of rare sources (N06, M08, lMiddelberg et alj|201lh . The AT- 
LAS survey regions were chosen because of the large amount of 
multiwavelength data covering these two fields, including near and 
far-infrared, and deep optical data, and in some areas, X-ray and 
UV. We aim to create the most comprehensive multiwavelength 
survey of faint radio sources to date. ATLAS currently contains 
~ 2000 radio sources, and we estimate we will have ~ 16000 ra- 
dio sources following final analysis of new observations completed 
in 2010 fro m ATCA with the Compact Array Broadband Back- 
end (CABB; IWilsonetalj201 1). These new observations consist of 
1000 hours of integration time, and cover a bandwidth of 0.5 GHz 
centred on 1.4 GHz. 

In ELAIS-S1, 1276 radio sources have been catalogued, 
comp rising 1366 radio components, at 1.4 GHz ( Mid delberg et al.l 
2008, M08). The 2.3 GHz image has a lower sensitivity, with an 
rms noise level of ~ 60/iJy in the central ~ ldeg 2 and ~ 100/iIy 
over the entire field. We find only 576 radio sources have a well- 
defined counterpart at 2.3 GHz, mainly because the resolution is 
~ 3 times coarser at this frequency than at 1.4 GHz as a conse- 
quence of the compact configuration of ATCA used for these ob- 
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Figure 2. (a) Grayscale MOST 843 MHz image of ATELAIS J003306.30-43 1029.8, (b) original ATCA 1.4 GHz greyscale image, and (c) the original ATCA 
2.3 GHz greyscale image. The white cross indicates the 1.4 GHz radio position of this source. 



servations. For the pu rposes of this pape r and catalogue, we use 
the cross-matches bv lZinn et alj Jin prepl Zl 1), where we take the 
1.4 GHz sources matched to the poorer resolution 2.3 GHz data, 
and extract the fluxes from these images for each source. M08 also 
compared their flux densities and positions to G99 by repeating 
their source extraction on the G99 1.4 GHz radio image. The dif- 
ferences in radio positions were determined to be negligible, but 
the flux densities of M08, while within ~ 3% of the earlier mea- 
surements, were consistently higher than G99. 



3 OBSERVATIONS, DATA REDUCTION AND 
CROSS-MATCHING 

To obtain low-radio frequency (843 MHz) data within ATLAS, 
we used the Molonglo Obs ervatory Synthesis Telescope (MOST; 
iMillsl 198 ^Robertson!! 199 ll) . The MOST is a 1.6 km long cylindri- 
cal paraboloid reflector, with its axis aligned in an east-west direc- 
tion. It operates at a freque ncy of 843 MHz, and has a field of view 
of 2?7cosec(<5)x2?7deg 2 jLarge et alJl98ll ; lBock et aljl999h . We 
have 31 separate 12 hour observations taken with MOST, which 
were combined into a single image. The CDFS field (centred at 
RA = 03 ,l 32 m 28 1 ' and Dec=-27°48'3l'.'0 (J2000.0)) was not ob- 
served, due to strong radio frequency interference (RFI) that in- 
creased with lower elevation. 



3.1 Processing 

The data reduction pipeline for MOST data is highly automated 
and reliable. Unfortunately, the telescope was subject to severe RFI 
at the time of our observations. The process to remove the RFI re- 
quired manual identification and excision of the affected data. Once 
completed, we were able to re-run the data pipeline, and produce 
final cleaned images. Seven of our 31 observations were subject to 
RFI that we were unable to remove completely, and these were not 
included in our final imaging. The remaining 24 images were added 
together in MIRIAD, first by re-gridding the images to ensure a 
common astrometry and pixel grid, and then using the task IM- 
COMB to combine the images. The final image has an rms sensitiv- 
ity of « 0.6 mjy, a higher value than expected due to the presence 
of low-level RFI that could not be excised. We catalogue sources 
down to a 5cr level of 3 mJy. The final image is shown in Figure Q] 
with the approximate borders of the ATLAS 1.4 and 2.3 GHz mo- 
saics overlaid. 

The MOST image was subject to two kinds of artifacts, grating 
rings and radial spokes. Due to the nature of the telescope design 



and observing mode, these artifacts are technically difficult to re- 
move. The artifacts remain in our image, and due care was taken to 
ensure all objects in our final catalogue were not spurious sources 
associated with the artifacts. An in -depth discussion of these arti- 
facts is given in lMauch et alj J2003I) . 

3.2 Source-finding 

Source finding for the final M OST image was done using the task 
SFIND jHopkins et alj|2002l) in MIRIAD. This produces a cata- 
logue of source positions, peak and total flux density and errors, 
plus the attributes of the Gaussian fits for each source. The final 
catalogue contains 325 radio sources, after removal of ~ 100 spu- 
rious sources associated with the grating rings and radial spokes. 
The spurious sources were identified and removed by visual inspec- 
tion (most lay directly on a strong radial spoke or diffraction ring), 
and by co mparison to the Sydney University M olonglo Sky Survey 
(SUMSS; lBock et al.ll999MMauch et alj2003h images (see ffHIfor 
details). Each MOST ATLAS source was identified in the SUMSS 
images with a low signal-to-noise source that was below the detec- 
tion level of the SUMSS catalogue. The MOST image is larger than 
the ATLAS survey region, so our catalogue encompasses a larger 
number of sources in total than we subsequently analyse together 
with the other ATLAS data. 

We use the process outlined in lHopkinsetai] ( l2003l) to calcu- 
late the errors associated with our m easured flux densities, and use 
Equation 5 of lHopkins et al. I d2003h with one small modification: 



07 =lJ 2.5^ + 0.05 2 , (1) 

where I is the total integrated flux density, a is the rms error in 
the image at the source location, and 07 is the total error on the 
integrated flux density. We use 0.05 2 as the sum of the squares of 
the instrumental a nd pointing er r ors a s given in M08, instead of 
0.01 as given in Hopki ns et al. I J2003h . We have used the more 
conservative estimate of these errors as given by M08, due to the 
calibration accuracy of MOST being comparable to that of ATCA. 

3.3 Cross-matching to ATLAS 

The final MOST catalogue was positionally cross-matched to the 
combined 1.4 and 2.3 GHz catalogue produced by Zll, that in- 
cluded the relevant SWIRE and optical data. As mentioned pre- 
viously in m.2\ only 576 of the 1276 1.4 GHz radio sources in 
ELAIS have a reliable single 2.3 GHz counterpart, whilst 460 have 
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Figure 3. (a) The greyscale MOST 843 MHz image of ATELAIS J003306.30-431029.8, (b) the convolved ATCA 1.4GHz greyscale image, and (c) the 
convolved ATCA 2.3 GHz greyscale image. The white cross indicates the 1.4 GHz radio position of this source, and the 1.4 GHz contours are overlaid at levels 
of 2.5, 5 and 12 mjy 



no counterpart, and 240 are blended or poorly fitted sources (where 
the 2.3 GHz counterpart encompasses multiple 1.4 GHz sources). 
From our cross-matching, 105 MOST sources were matched to sin- 
gle sources in the Zll catalogue. Another 30 MOST sources were 
matched to confused or blended sources from Zll, and 31 MOST 
sources were found to have multiple isolated Zll sources matched 
to a single MOST source. 



3.4 Resolution Matching 

The resolution of the MOST image, 62" x 43", is much coarser 
than the resolution of the 1.4 GHz image (10" x 7"), and the 
2.3 GHz image (33" x 20"). To determine accurate spectral in- 
dices across the three frequencies, it is necessary to convolve the 
ATLAS 1.4 and 2.3 GHz images to the same size as the MOST 
beam to ensure the recovered flux for sources at all frequencies in- 
cludes any emission extended on scales up to those consistent with 
the MOST beam. Figure [2] shows an example of a complex source 
from the ATLAS 1 .4 GHz catalogue in greyscale in the three radio 
frequencies, to highlight the differences in resolution. In Figure [3] 
the same source is shown after convolving the 1.4 and 2.3 GHz im- 
ages to the same resolution as the MOST image. Although convolu- 
tion removes most of the small scale structure visible in the original 
ATLAS 1.4 GHz image, it ensures we are detecting emission from 
the same spatial region from each source. 

For consistency, we used SFIND to detect sources within the 
convolved 1.4 and 2.3 GHz images. This allows us to measure er- 
rors on our flux densities that are consistent for the three images, 
which is important for producing accurate spectral indices and as- 
sociated uncertainties. Rather than rejecting blended sources en- 
tirely, we account for them carefully in our catalogue. These are 
single or point sources in the 843 MHz MOST image (and the 1.4 
and 2.3 GHz convolved ATLAS images) that encompass multiple 
ATLAS 1.4 GHz sources from the M08 catalogue (Figure g}. The 
M08 flux densities for these sources reported in our catalogue are 
the summed values of the individual M08 1.4 GHz flux densities. 
As the contribution to the total flux density from each individual 
component within a blended source cannot be determined, we treat 
them as a single entity for the purpose of the current analysis, even 
though they may be physically unrelated. We include these sources 
in our analysis for completeness, but caution that their spectral in- 
dices should not be considered as accurate estimates of those for 
the underlying source components. 



3.4.1 Flux Density Comparisons and Corrections 

To ensure our measured 1.4 and 2.3 GHz flux densities from the 
convolved, coarser-resolution, images were robust, we compared 
these to the flux densities in M08 and Zl 1. We find that our mea- 
surements in the images convolved to a resolution consistent with 
the MOST image tend to systematically overestimate the 1.4 GHz 
flux densities compared to M08, and underestimate the 2.3 GHz 
flux densities compared to Zl 1. 

The issue at 2.3 GHz arises from negative CLEAN bowls 
around sources in the original image. These artifacts give rise, after 
convolution, to our observed systematic decrease in flux density. 
The 2.3 GHz beam size (33" x 20") is close to that of our MOST 
images (62" x 43"), and the bulk of the sources are unresolved, 
meaning that there is likely to be little flux missed on extended 
scales. Consequently, we choose to use the flux densities estimated 
by Zll in the original 2.3 GHz image for our spectral index esti- 
mates. 

At 1.4 GHz, the difference in resolution from the MOST im- 
age is sufficient that we need to use the flux densities from the 
convolved image. The flux density overestimate compared to M08, 
which is limited to the fainter sources (5843 < 10 mjy), is typi- 
cally of the order of 10%. This exists for clearly unresolved sources, 
which should be identical before and after convolution. This issue 
is associated with side-lobes from the radio sources in each individ- 
ual telescope pointing, below the level to which the image has been 
CLEANed, being summed in the convolved image. While this is a 
small effect, with a minimal impact on our derived spectral indices 
(quantified below), we can make an empirical correction that mini- 
mizes any impact further still. This is implemented through a least- 
squares fit of our convolved flux densities against those of M08, and 
scaling our convolved flux densities using this fit to be consistent 
with those of M08. This accounts for the flux-density dependence 
of the imaging systematics while retaining contributions from any 
real extended flux components, and at the same time, minimizing 
any possible overestimate of the flux density. This correction typi- 
cally results in a change of only ctf it ~ 0.02. 

The impact on our derived two-point spectral index estimates 
of potential remaining flux density uncertainties at 1.4 and 2.3 GHz 
of ~ 10% is ~0.1. This uncertainty is included in our estimates of 
spectral index errors below, by being added in quadrature with the 
flux density errors (given by Equation QJ. 
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Figure 4. An example of a blended source. The contours are from the con- 
volved 1.4 GHz image (at levels of 5, 10 and 20mJy), and the greyscale is 
the 843 MHz image. The white crosses indicate the 1.4 GHz radio position 
of the two sources, where the object in the centre of the image is the MOST 
cross-matched source, and the blended object indicated by the other cross. 



3.4.2 Blended sources 

Another effect of the convolution is that most objects classified as 
radio doubles, triples, or core-jet morphology by M08, appear as 
one MOST object (single or blended). Only two core-jet morphol- 
ogy sources appear slightly elongated in the direction of the jet in 
the MOST and convolved 1.4 and 2.3 GHz images. 

Due to the lower resolution of the MOST image, there are 61 
instances where there are ATLAS sources within the beam for a 
MOST radio source. These 61 MOST sources correspond to a to- 
tal of 144 ATLAS M08 1.4 GHz sources. Of the 61, 30 sources 
were classified as blended by Zll due to having multiple 1.4 GHz 
sources within the 2.3 GHz beam (33" x 20"). The remaining 31 
are classified as blended in this paper due to multiple 1 .4 or 2.3 GHz 
sources within the MOST beam. The cross-matched 1.4 GHz AT- 
LAS source (positionally closest to the MOST radio position) is 
generally bright while the second, blended, source contributing to 
the flux density is typically much fainter (Sia < 1 mjy). As a con- 
sequence, the majority of such blended sources are likely to have 
little contribution to the MOST flux density from the secondary 
component, and the spectral index estimate for the blended sources 
can be considered as that for the primary ATLAS counterpart, even 
though the uncertainties on this estimate will clearly be larger than 
the formal uncertainties provided in our catalogue. Blended sources 
are flagged in the catalogue. 

For blended sources in our convolved images, which were not 
classified as blended sources by Zl 1, the flux densities from sepa- 
rate sources in the ATLAS Zll catalogue are summed. A source is 
defined as being part of a blended object if the 1.4 GHz radio posi- 
tion lay within an ellipse the size of the Gaussian used to find the 
total flux density in the convolved image. An example of a blended 
source is shown in Figure|4]where the cross-matched ATLAS coun- 
terpart is in the centre of the MOST source, and the blended object 
is within the MOST synthesized beam. There is one exception to 
this, shown in Figure [5] where the blended object is outside the 
Gaussian fit, but the radio emission clearly extends from S120 into 
S107. 



O -23' 




Yarns' 



60" 45" 

Right Ascension (J2000) 



Figure 5. The single example of a blended source where the ATLAS 
1 .4 GHz radio position is not within the Gaussian (shown by the white el- 
lipse) fitted to this source, but there is clearly contamination from the neigh- 
bouring source. 



4 THE MOST ATLAS SOURCE CATALOGUE 

The total catalogue consists of 325 MOST sources, limited to 
sources above the 5a cutoff of 3 mjy, of which 166 have an ATLAS 
1.4 or 2.3 GHz counterpart. The remaining 159 sources are outside 
the ATLAS survey area. In total, there are 310 ATLAS sources in 
our catalogue, with 205 classified as blended sources. The remain- 
ing 105 ATLAS sources correspond to single MOST sources. An 
extract of the catalogue in shown in Table [TJ and the full version 
is available online. The MOST radio position is listed first, fol- 
lowed by the flux density and error at 843 MHz. For MOST sources 
with a single ATLAS counterpart, the ATLAS source name (e.g. 
ATELAIS J002905.22-433403.9), and ID (e.g. S100) is listed; 
for blended sources the given ATLAS IDs correspond to all Zll 
sources cross-matched to the MOST source, and the ATLAS source 
name corresponds to the first ATLAS ID listed. All ATLAS source 
names and IDs are from M08. The corrected convolved 1.4 GHz 
and Zll 2.3 GHz flux density measurements and associated errors 
are given, along with the three spectral indices (described in 
below). 



5 RESULTS AND ANALYSIS 
5.1 Flux Density Distribution 

The flux density distributions (Figure [6]l do not show any major 
differences between the single, blended, and non-ATLAS sources. 
The location of the median of the entire catalogue is shown in 
Figure |6jt. The distributions appear consistent with ot her faint ra- 
dio samples, such as the orig inal ATLAS catalogues ( Norr is et al.l 
l2006tlMiddeiberg et alj2008h . 

The Sydney University Molo nglo Sky Survey (SUMSS; 
iBock et alj[l999l iMauch et alj|2003h is an 843 MHz survey with 
MOST, that covers the sky south of 6 < -30° with \b\ > 10°. 
SUMSS has similar resolution and sensitivity to the National Ra- 
dio Astronomy Ob servatories (N RAO) Very Large Array (VLA) 
Sky Survey (NVSS: ICondon et al.ll998h . SUMSS contains 21 1063 
radio sources, with an rms of ~ 1 mJy. Our observations probe 
sources fainter than SUMSS by a factor of ~ 2. Of our catalogue, 
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Figure 6. (a) The flux density distribution for our entire catalogue of MOST sources with the median of the catalogue as the vertical dashed line, (b) Distribution 
of flux densities for single and blended MOST ATLAS sources, and (c) the MOST sources outside the observed ATLAS ELAIS field. 
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Figure 7. (a) Spectral index a f it versus 843 MHz flux density, (b) 1 .4 GHz flux density, and (c) 2.3 GHz flux density. The dotted black line indicates the ltr 
spectral index limit, and the solid black line represents the 5cr spectral index limit for each flux density. The dashed black line is the median a fa value for 
single sources only. 
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Figure 8. Distribution of a f it for our catalogue. 



178 sources are cross-matched to SUMSS sources. Only one source 
in SUMSS does not have a counterpart in our observations, as it is 
located directly on an artifact in our image. Our measured flux den- 
sities, positions, and radio differential source counts (see £15.51 for 
details) are consistent with those of SUMSS, with only two out- 



liers, likely due to intrinsic source variability. It is known that a few 
percent of mjy radio sour ces are variable on the timescale of years 
dOort & Windhorstlll985h . Attributing our outliers to variability is 
consistent with this work, as our data was taken over several years. 



5.2 Spectral Index Distributions and Properties 

We have calculated spectral indices for all sources for which flux 
densities are available at two or three frequencies, a f it is a three- 
point power-law fit, and Qo!843 an d a i'.4 are two-point power-law 
fits for the respective frequencies. Distributions of these spectral 
indices are given in Figures |7j and [8] Spectral index ao!s43 versus 
spectral index af;| is also shown in Figure [9] which indicates that 
most of our sources are steep-spectrum, with a small fraction of 
inverted, peaked, and flat spectrum objects, discussed further be- 
low. A large proportion of steep-spectrum sources are seen because 
of the low radio frequency selection, and steep spectrum sources 
are brighter at lower frequencies. The median spectral index a fa 
for the single sources is shown on each of the panels in Figure [7] 
Only single sources are included in this calculation. The spectral 
index limits are shown in Figure|7p,c arising from the limiting flux 
density of the least-sensitive frequency in the relevant two-point 
spectral index calculation. Although the three-point spectral index 
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Table 1 . Extract from the 843 MHz MOST ATLAS Catalogue 



oo 



MOST RA 


MOST Dec 


So. 843 


A So. 843 


ATLAS Name 


ATLAS ID/s 


Type 


Si. 4 


ASi.4 


S2.3 


AS 2 . 3 




A a fit 


1.4 
0.843 


1.4 

a 1 3 
2.3 






mJy 


mJy 








mJy 


mJy 


mJy 


mJy 




0:28:45.438 


-42:51:39.46 


15.12 


1.5 
























0:28:54.063 


-43:12:18.30 


6.80 


1.8 
























0:28:56.980 


-42:18:15.37 


40.52 


3.2 
























0:29:05.054 


-43:34:07.05 


11.34 


1.3 


ATELAIS J002905.22-433403.9 


S749 


P 


9.01 


1.16 


5.1 1 


0.34 


-0.81 


0.06 


-0.45 


-1.14 


0:29:05.968 


-44:42:00.59 


14.18 


1.6 
























0:29:09.167 


-43:44:02.10 


11.13 


1.3 


ATELAIS J002909.26-434356.3 


S617 


P 


12.52 


1.27 


4.63 


0.27 


-0.93 


0.09 


0.23 


-2.00 


0:29:13.163 


-44:52:23.32 


5.33 


1.4 
























0:29:15.437 


-43:26:36.78 


9.51 


1.3 


ATELAIS J002915.52-432638.3 


S868 


P 


5.45 


1.10 


3.56 


0.25 


-0.97 


0.03 


-1.10 


-0.86 


0:29:21.471 


-42:55:45.26 


44.87 


2.6 
























0:29:25.667 


-44:08:25.68 


15.37 


1.8 


ATELAIS J002925.66-440822.8 


S293, S304 


b 


9.12 


1.13 


3.51 


0.35 


-1.49 


0.06 


-1.03 


-1.92 


0:29:27.845 


-43:16:15.59 


6.61 


1.2 


ATELAIS J002927.69-431614.4 


S1014 


P 


3.85 


1.16 


2.03 


0.18 


-1.18 


0.01 


-1.06 


-1.29 


0:29:36.620 


-42:25:41.20 


64.07 


8.4 
























0:29:37.206 


-42:34:18.69 


5.05 


1.4 
























0:29:38.084 


-44:23:21.78 


42.38 


2.7 


ATELAIS J002939.19-442319.3 


S100, S101.1 


b 


34.1 


2.29 


18.82 


1.88 


-0.74 


0.06 


-0.43 


-1.20 


0:29:43.820 


-42:37:47.77 


17.83 


1.7 
























0:29:45.687 


-43:21:50.03 


24.18 


1.6 


ATELAIS J002945.64-432149.5 


S943 


P 


16.25 


1.41 


10 


0.52 


-0.87 


0.06 


-0.78 


-0.98 


0:29:46.493 


-43:15:57.34 


41.67 


2.4 


ATELAIS J002946.52-431554.5 


S1018 


P 


27.89 


1.78 


18.49 


0.94 


-0.81 


0.05 


-0.79 


-0.83 


0:29:47.622 


-44:16:15.59 


6.62 


I.I 


ATELAIS J002947.37-441607.0 


S181 


P 


3.90 


1.65 


1.75 


0.15 


-1.33 


0.01 


-1.04 


-1.62 


0:29:50.398 


-44:05:48.45 


7.16 


I.I 


ATELAIS J002949.89-440541.4 


S345, S342, S339 


b 


5.26 


1.07 


5.13 


0.5 1 


-0.31 


0.03 


-0.61 


-0.05 


0:29:51.452 


-43:45:28.25 


8.23 


1.3 


ATELAIS J002951.48-434528.0 


S598 


P 


4.93 


1.13 


3.23 


0.18 


-0.93 


0.02 


-1.01 


-0.85 


0:29:53.062 


-42:50:17.87 


5.16 


2.0 
























0:30:03.253 


-42:18:22.97 


20.81 


8.4 
























0:30:04.749 


-42:09:58.88 


19.10 


5.9 
























0:30:10.595 


-44:09:12.06 


6.49 


I.I 


ATELAIS J003010.82-440907.3 


S288 


P 


7.60 


0.72 


8.96 


0.46 


0.32 


0.04 


0.31 


0.33 


0:30:17.439 


-42:24:47.35 


436.20 


22.4 
























0:30:18.530 


-45:26:42.05 


12.79 


2J 
























0:30:18.793 


-44:04:35.34 


13.70 


1.2 


ATELAIS J003019.22-440438.3 


S355 


P 


9.83 


0.78 


5.57 


0.23 


-0.89 


0.0s 


-0.65 


-1.14 


0:30:20.948 


-43:39:44.51 


69.27 


3.7 


ATELAIS J003020.95-433942.8 


S694 


P 


51.99 


2.55 


33.92 


1.7 


-0.66 


0.04 


-0.57 


-0.86 


0:30:21.663 


-45:05:09.55 


16.15 


3.3 
























0:30:22.709 


-42:37:02.29 


14.25 


2.4 
























0:30:27.017 


-42:35:14.14 


7.34 


3.1 
























0:30:29.096 


-42:13:50.96 


18.19 


6.2 
























U. JU.Jt.ojj 


— HO.ZV.H- / .L\> 




























0:30:35.339 


-44:37:11.21 


3.69 


I.I 


ATELAIS J003035.77-443707.2 


SI8.S19 


b 


6.53 


3.21 


2.07 


0.21 


-0.68 


0.01 


1.13 


-2.31 


0:30:35.937 


-43:23:39.75 


8.95 


3.1 


ATELAIS J003035.03-432341.6 


S926, S923, S930, S930.1 


b 


5.03 


0.62 


4.76 


0.48 


-0.32 


0.02 


-1.14 


-0.11 


0:30:38.957 


-44:09:56.34 


3.95 


0.9 


ATELAIS J003039.03-441000.0 


S279 


P 


3.04 


0.62 


1.23 


0.13 


-1.29 


0.01 


-0.52 


-1.82 


0:30:39.015 


-45:07:07.04 


33.29 


2.7 
























0:30:39.621 


-44:42:01.21 


28.43 


1.9 


ATELAIS J003039.68-444159.5 


S7 


P 


26.00 


3.44 


17.02 


0.95 


-0.49 


0.05 


-0.18 


-0.85 


0:30:40.860 


-43:23:40.55 


11.17 


2.4 


ATELAIS J003042.10-432335.4 


S923, S930, S930.1 


b 


6.02 


0.64 


5.47 


0.55 


-0.53 


0.03 


-1.22 


-0.19 


0:30:42.041 


-43:18:42.39 


4.52 


1.0 


ATELAIS J00304I.88-431840.7 


S987 


P 


3.26 


0.61 


2.12 


0.13 


-0.77 


0.01 


-0.64 


-0.87 



y-jit is the spectral index fitted across the three flux density measurements at 0.843, 1.4 and 2.3 GHz. Type (Column 7) refers to whether the source is a single point source (p) at all frequencies, or a blended source (b). 
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Table 2. Spectral Classifications 



Type 


Total 


Point 


Blended 


Percentage 




Number 


Sources 


Sources 




Steep-Spectrum 


143 


89 


54 


86% 


Inverted 


3 


2 


1 


2% 


Peaked 


12 


8 


4 


8% 


Upturn 


6 


4 


2 


4% 



NOTES.- The percentages of each source are only for the single or point 
sources. 



am is shown in the figure, the limit associated with the two- 
point spectral index calculation provides a clear indication of where 
we are selection-limited against particularly steep or flat spectrum 
sources. Distributions of ao!843> split into three flux density bins 
are als o shown in Figure [10] analogous to Figure 8 of lMauch et al.l 
(2003). A flattening of the spectral index «q;| 43 is suggested by 
Figure ITol however this is primarily due to different numbers of 
sources in each flux density bin. This is further discussed in ^5 .41 
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Figure 9. Spectral index Oj | versus for MOST ATLAS sources. 

The black symbols in each outer corner of the plot represent the type of 
source in each quadrant of the plot, e.g. steep, inverted, peaked or upturn. 



5.3 Radio spectral classifications 

Figure|9]sh ows our sample of sources, in a radio colour-colo ur dia- 
gram jKesteven. Bridle & Brandidl 19771 : llVlurphv et alj|2010h . that 
we now consider in four classes: 

(i) Steep-spectrum objects, with a steep radio spectrum from 
843 MHz to 2.3 GHz (the lower-left quadrant of Figure[9}, 

(ii) Peaked sources, where we see the radio spectrum turn over 
between 843 MHz and 2.3 GHz (the lower-right quadrant of Fig- 
ured],, 

(iii) Inverted (or rising) sources, where the radio flux density 
increases with increasing frequency (the upper-right quadrant of 
Figure[9],, and, 

(iv) Upturn sources, that have an upturn in their radio spectrum 
(the upper-left quadrant of Figure[9],. 

This distribution highlights the fact that our sample is dominated 
by steep-spectrum sources. The statistics are given in Table [2] 
split into single and blended sources, noting that any statistics 
from the blended sources are not reliable. In comparison, the Aus - 
tralia Telescope 20 GHz Survey (AT20G; iMurphv et af] l2010h . 
a high-frequency- selected sample (20 GHz), produced a colour- 
colour plot contained 3763 sources, with 14% inverted (rising) 
spectrum sources, 57% steep-spectrum objects, 21% peaked spec- 
trum sources, and 8% sources with an upturn in their spectra. 

As indicated by Tabled 8% of our MOST ATLAS sample 
appear to be possible Giga hertz Peaked S pectrum (GPS) sources 
(see Sj6T]for a description). lO'Deal ( fl998h su ggest that ~ 10% of 



bright radio sources are GPS sources, whereas Randall et alj feOllh 
found less than 1% of their sample to be GPS sources. Although 
both these samples are not complete, it is interesting to note that at 
faint radio fluxes, the propo rtion of GPS s ources appears to be sim- 
ilar to the bright sample of O'Deal J 19981) rather than the unbiased 
bright sample o f lRandall et alJfeOllh. This sample of sources will 
be explored further in a future paper, Randal fet al.l Jin prepj) . 

Interestingly, there is one very steep spectrum single source, 
S1256, which has «o!843 = —2.92, that is not detected in the 
2.3 GHz image. This object is discussed in more detail in i]7.3.5l 



Single Source: 
Blended Sourc 



(a) .S' ;i . J:1 <10mJy 



(b) 10mJy<S„ M:1 <20mJy 



-2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 
Spectral Index 

Figure 10. Distribution of aj | 43 , in three flux density bins, 0-10mJy (a), 
10-20 mJy (b), and > 20mjy (c). The solid vertical lines are the median 
spectral index ctj | 43 for each flux density bin for the single sources only. 



5.4 Spectral index as a function of flux density 

We have investigated the properties of our sample with flux density. 
FigurefTTIshows median spectral index as a function of flux density, 
both for our sample (Figure [TTb). and for a compilation of samples 
from the literature along with ours (Figure [TTb). While the uncer- 
tainties are large, there is evidence for a mild trend toward a flat- 
ter spectral indices with decreasing flux density (Figure [TTV). Fig- 
ure II lb shows the com p arison of our m e dian s p ectral indices wit h 
Wind horst et alj dl993l) : iPrandoni et ail ( l201lh : llbar et al.l J2009I) . 
and lOwen et al. (2009). Our data is in general consistent with these 
previous results, and so we cannot rule out the possibility that there 
is a flattening of the median spectral index with decreasing flux 
density. We note that the very steep median spectral index in the 
faintest flux density bin of many surveys is a consequence of the 
flux density limits preventing the detection of flatter spectral in- 
dices close to the survey limits. 
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Figure 11. (a) Median spectral index aj'| 4 3 for our sample with 843 MHz flux density, (b) Median spectral index '| 43 for our sample, compared to median 
spectral indices from the literature. The Windhorst et al. 1993 solid line indicates a fit to the median spectral index g for their data. The other values 

are taken from tabulated results in Ibar et al., 2009, Owen et al., 2009, Prandoni et al., 2011, and Mauch et al. 2011 (Priv. Comm). The horizontal error bars 
represent the flux density bin for each median point; the vertical error bars are the 25" 1 and 15 th percentiles. 



Table 3. Median Spectral Index ctj g43 Statistics as a function of flux den- 
sity 



Range in 5843 
or Si. 4 (mjy) 


Median 

"0.843 


Mean Flux 
(mjy) 


Number 
of sources 


25 th 
perc. 


75 th 
perc. 


Si. 4 


0.81-3.7 
3.7-6.4 
6.4-9.8 
13.5-26 
26-95.4 


-1.134 
-0.717 
-0.533 
-0.584 
-0.775 


2.50 
4.89 
7.68 
13.65 
36.21 


21 
21 

22 
21 
20 


-1.390 
-1.016 
-0.876 
-0.769 
-0.945 


-0.676 
-0.355 
-0.275 
-0.332 
-0.656 


•So. 843 


3.3-4.52 
4.52-6.7 
6.7-11 
11-15 
15-35 
35-204 


-0.818 
-0.488 
-0.676 
-0.671 
-0.687 
-0.921 


4.04 
5.54 
8.76 
12.68 
22.71 
62.42 


18 
18 
17 
18 
18 
16 


-1.273 
-1.203 
-1.210 
-1.045 
-0.872 
-0.930 


-0.825 
-0.733 
-0.817 
-0.795 
-0.575 
-0.741 



NOTES.- Errors on the median spectral indices are calculated from the 
25th and 75th percentile. The mean flux densities are calculated as the 
mean of the flux densities of all the sources in each bin. 



5.5 Radio Source Counts 

The differential radio source counts for our 843 MHz data are pre- 
sented in Table|4] A weighting factor has been applied to the source 
counts to correct for incomplet eness due to the nois e level increas- 
ing at the edges of the field (see Hopki ns"et alj 19981) . We do not ap- 
ply a resolution correction, as the MOST has a large be am, and has 
a high sensitivity to ext ended diffuse radio emission dBock et all 
1 1999i;lMauchetal .120031) . The differential source counts are shown 
in Figure [12] with the SUMSS jMauch et al] 120031) differential 
source counts, and a comp ilation of 1.4 GHz differential source 
counts (Hopkin s et al. 120031) as a reference sample. Our data probes 
the 843 MHz source counts a factor of 2 fainter than SUMSS, but 
we still clearly underestimate the counts in our faintest flux density 
bin, due to incompleteness. 

The reference differe ntial source count sample at 1.4 GHz 
from iHopkins et~al1 d2003l) has been shifted to 843 MHz assum- 
ing an average spectral index value. We have determined that 
a = —0.5 is the appropriate value to use for the average spectral 



index, using the Kellermann correction, discussed in further detail 
here. 



5.5.1 The Kellermann Correction 

iKellermanrj i ll 9641) noted that any observed spectral index distri- 
bution is not independent of the observing frequency. The correc- 
tion accounts for the relationship between observed spectral index 
distributions at different frequencies. At higher radio frequencies, 
we tend to see more flat spectrum objects as proportionally more 
sources are above the flux limit. In contrast, as we move to lower 
radio frequencies, we tend to see more steep spectrum objects be- 
cause they are brighter at the lower radio freque ncies. The correc- 
tion stated in the Appendix o f lKellermanrJ l l 19641) describes the off- 
set between mean values of spectral indices at different frequencies. 
We use Equation A5 of lKellermannl (1964), where we assume we 
have two distributions of spectral indices {P(a) and Q(a)), at two 
different frequencies, v\ and V2 respectively. Also requir ed is x, 
the sl ope of the number-flux power-law N(> S) = k S x iLongair] 
1966), where k is a constant scaling factor, and TV represents the 
number of sources above a given flux density, S. 

Q(a) = A(^)- a *P(a), (2) 

/oo 
Q(a)da = 1. (3) 
-oo 

P(a) is the known distribution of spectral indices at one frequency 
and Q(a) is the distribution of s pectral indices we wish to in- 
fer, given by Equation|2] Following lKellermannl d 19641) we assume 
the special case where these two distributions are Gaussian (a rea- 
sonable assumption given the shapes of the distributions in Fig- 
ures [8] and [TUJ, and have t he same dispersion a, and A is chosen 
such that Equation A6 of lKellermannl i 19641) (Equation [3} is sat- 
isfied. This special case scenario results in the mean value of the 
inferred spectral index distribution (Q(a)) being shifted by a fac- 
tor of xo 2 h\(vi j i*2 ) • For our sample, to shift the distribution of 
the 1.4 GHz source count compilation to 843 MHz, we used aj;| 43 
to determine the value of the Kellermann correction for only the 
single sources. The dispersion was found to be a = 0.56 and 
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Table 4. Source Counts from the ATLAS ELAIS MOST Observations 



Range in 5g43 


N 


Neff 


( "5843 ) 


(dN/dS)/S~ 2 - 5 


(mJy) 






(mJy) 


(Jy^sr" 1 ) 


3.31-5.81 


45 


21.58 


4.38 


4.29 


5.81-8.31 


45 


34.87 


6.95 


21.90 


8.31-10.81 44 


41.63 


9.48 


56.83 




10.81-15.81 


51 


50.02 


13.07 


76.29 


15.81-25.81 


48 


47.54 


20.20 


107.61 


25.81-45.81 


42 


42 


34.38 


179.70 


45.81-100 


33 


33 


67.68 


283.27 



N e f f is the effective number of sources in each bin, after the weighting 
factor has been applied. The mean flux density given is the geometric 
mean of the bin limits. 



-a ict 



□ 1.4 GHz Source Counts Compilation (Hopkins et al. 2003) shifted to 843 MHz by c 
4 SUMSS Source Counts (Mauch et al. 2006) 
• MOST ATLAS (this paper) 




1CV 1 
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Figure 12. The differential source counts for our data (circles), SUMSS (di- 
amonds), and a 1 .4 GHz compilation of source counts (grey squares) shifted 
to 843 MHz assuming a spectral index of a = —0.5, from Hopkins et al. 
2003 as our reference sample. 



x = —1.60, resulting in a Kellermann correction of 0.25. This 
shift in the mean of the spectral index distribution from 843 MHz 
to 1 .4 GHz is enough to account for the use of the a ~ —0.5 to shift 
the 1.4 GHz source counts to the 843 MHz source counts, instead 
of the observed mean spectral index of our sample a — —0.71. 



6 FAINT CSS AND GPS CANDIDATES SOURCES IN 
ATLAS 

6.1 Compact Steep Spectrum and Gigahertz Peaked 
Spectrum Sources 

Gigahertz Peaked Spe ctrum (GPS) and Compac t Steep Spec- 
trum (CSS) sources dPolatidis & Conwavl 120031 ; iFantil 12009a; 
Morg anti et ail ]2009) are a class of compact, powerful radio 
sources, suggested to be the beginning of the evolutionary path for 
large-scale radio sources. The proposed evolutionary path has GPS 
sources evolving into CSS sources, whic h then gradually evolv e 
into Fanaroff-Riley Type I and II galaxies dFanaroff & Rilevl 1974b . 
depending on their initial luminosity. These sources offer an ideal 
resource to investigate galaxy evolution and formation, as well as 
AGN feedback, as they are young AGN, but also have star for- 



mation occurring due to interactions an d mergers dO'Dea|[l99ot 
Labia no et al . 20081 iMorgantT et al. 2009). A mor e detail ed discus- 
sion of their properties is given in Randall et al. U2011I). If bright 
CSS sources do evolve into FRI/IFs dFanaroff & Rilevll 1974b. then 
some strong evolution must occur for this to happen dSnellen et alj 
1999), as we do not see this high percentage of CSS and GPS 
sources in the local Universe. It is also possible that the supply 
mechanism of the energy powering these objects could cut off, 
leaving only diffuse emission and a steep spectrum core. This 
would result in faint objects whose radio l obes have c eased to ex- 
pand or a prematurely dying radio source ( Fanti 2009b), explaining 
the lack of large numbers of these objects nearby. 

Bright CSS and GPS sources are common (~ 30% and 10% 
respectively) in radi o surveys, bu t few f a int samples exi s t. Previ - 
ous surveys inclu de Snelle n et alj d2000]) ; lTschager et alj { 2003b); 
iKunert-Bairaszewska & Mareckil d2007l) ; iFantietalJ d201ll and 
references therein;), which catalogue CSS and GPS sources down 
to ~ 20 mJy. If these objects are as prevalent at faint flux densi- 
ties, a whole population of these objects r emains mostly unknown 
dSnellen et al.ll998tlTschager et al. 2003a). Here we present an ini- 
tial complete sample of faint CSS and GPS sources from ATLAS, to 
further our understanding of their role in galaxy formation and evo- 
lution. Understanding their properties across different wavelength 
regimes is important, as is studying samples across a wide range 
of flux densities, as this will help to build a complete picture of 
their properties and nature. We present here a brief overview of the 
selection and properties of the faint sample of candidate CSS and 
GPS sources selected from ATLAS. 



6.2 A new candidate sample of faint CSS sources 

An initial sample of faint CSS sources has been selected from AT- 
LAS ( based upon spectral index information, and angular 
size. The initial sample is drawn from both ATLAS fields, CDFS 
and ELAIS, and only includes the 1.4 and 2.3 GHz flux density 
measurements. CDFS was not observed by MOST, and we do not 
consider the ELAIS MOST observations in this discussion. 

The selection criteria listed in Table[5]were utilized to first se- 
lect a sample of unresolved, single objects from both the ELAIS 
and CDFS fields, before applying a spectral index cut. The number 
of sources remaining in each field after each criterion is also listed. 
We have chosen a primary spectral index cut of af;| < —0.95 
with an additional supplementary set extending to af l < —0.9, 
as the errors on the spectral indices are typically ~ 0.2. These 
values were chosen to avoid selecting any star forming galaxies 
(with typical synchrotron spectral indices of af ^ ~ —0.7), as 
at these low flux dens ities, star forming galaxies are a dominant 
part of the population dWindhorst et alj|l985l 1 19931 ; iHopkins etall 
l2003l ; ISevmour et alj|2008b . In the ELAIS catalogue, there are 576 
sources with a single 2.3 GHz match (blended sources are not in- 
cluded in this analysis) out of 1276, and 399 of 726 in the CDFS. 
The primary reason for only half the sources having a 2.3 GHz 
counterpart is the lack of equivalently deep 2.3 GHz data across 
the two fields. The 2.3 GHz data has lower resolution (33" x 20" 
and 54" x 21" for ELAIS and CDFS respectively), and does not 
have the same sensitivity as the 1 .4 GHz data for both ELAIS and 
CDFS. The rms of the CDFS 2.3 GHz data is ~ 3 times higher than 
the ELAIS 2.3 GHz data. The initial selection criterion of being 
an unresolved component also removes 63% of the sources from 
CDFS, but only 47% from ELAIS. The percentage of sources re- 
maining after the initial selection is small, with 86 sources and 5 
supplementary sources, from a possible 975 (~ 10%). 



© 2011 RAS, MNRAS 000. [71011 



12 Randall et al. 

Table 5. Summary of selection criteria and number of sources remaining 



Field 




ELAIS 


CDFS 


Criterion 


Selection 


Sources 


Sources 


Number 


Criterion 


Remaining 


Remaining 


1 


Unresolved in component catalogue 


681 


271 


2 


Not part of a radio double/triple/complex object or classified as a sidelobe 


638 


247 



3 Apply spectral index cut of | < —0.95 for sample selection 78 8 

3a Sources with spectral index between —0.9 > > —0.95 selected as supplementary sample 4 1 

NOTES. -Removal of radio doubles/triples/complex sources was done by cross-matching both the component, and source catalogues from CDFS and ELAIS, 

as the source catalogue lists what components comprise each source. 





Figure 13. (a) The 1.4 GHz M08/N06 flux density distribution for the initial CSS sample, for the full and supplementary sets, shown as the grey filled 
distribution. The entire ATLAS catalogue is shown in the solid empty distribution. The solid vertical line represents the median of both distributions, (b) 
The solid grey filled histogram indicates the spectral index oJ,'! distribution for the initial CSS sample, and the solid empty distribution is the entire ATLAS 
catalogue spectral index distribution. The vertical solid and dashed lines indicate the median of each distribution respectively. 



Of the 86 sources, only 5 are cross-matched to 843 MHz 
sources, due to the relatively high flux limit of the 843 MHz data. 
The relevant information from these sources is added into the sam- 
ple catalogue, but is not used for any subsequent analysis below. 
In considering the spectral index selection, described in Table [5] 
for each field, we note that the CDFS 1 .4 GHz flux densities were 
taken for the sources as observed, rather than after convolution to 
match the 2.3 GHz resolution. This step has been omitted from the 
current analysis due to the small numbers of sources involved, and 
while the spectral index estimates for these few sources may be 
somewhat less robust than those in the ELAIS field, this does not 
have a significant impact on any of the results presented here. 

In the ELAIS field, we use the spectral indices measured be- 
tween 2.3 GHz and the Zll 1.4 GHz flux densities (measured with 
the same resolution as the 2.3 GHz image). If the source does not 
have a reliable convolved measured flux density, we use the M08 
flux (this is only necessary for four sources). To robustly constrain 
the radio spectra of these sources, we need to explore both higher, 
and lower frequency data to confirm that our sources are truly steep- 
spectrum or peaked spectrum objects. We show the basic statis- 
tics of our sample in Table|6] and discuss initial results in i|7] The 
1.4 GHz flux density and spectral index distributions are presented 
in Figure [T3] with the median flux density and spectral index in- 
cluded as the solid black lines. The spectral index distribution of 
the entire ATLAS catalogue is also shown in Figure [13b to provide 



context for the initial CSS sample distribution. Note that the steep 
spectrum sources not identified as CSS are those that were rejected 
by criteria 1 or 2 in Table [5] The median of the CSS sample is 
a\ : X = -1.48, compared to a\'\ = -0.84 for the entire ATLAS 
catalogue. 



7 DISCUSSION 

7.1 Redshift and Luminosity Properties 

ATLAS has ^700 spectroscopic dMao et alJIin prep.1) . and ~700 
photometric dRowan-Robinson et a D l2004h redshifts for radio 
sources in the CDFS and ELAIS-S 1 fields. The sources chosen for 
the spectroscopic targets were primarily the optically bright sources 
(80% complete at R magnitude ~20), that are mainly low redshift 
sources, although some quasars at higher redshift are also selected. 
For the 105 single sources in the MOST ATLAS catalogue, we 
found 20 spectroscopic redshifts and 9 photometric redshifts (no 
analysis is done on the blended sources). In the initial CSS sample 
of 91, 22 spectroscopic and 17 photometric redshifts were found. 
In Figure[l4]spectral index against redshift and 1.4 GHz rest-frame 
luminosity are shown, for both the MOST ATLAS catalogue and 
the initial CSS sample. The grey box in Figure [14b indicates where 
SFGs would most likely be located, defined by the dividing line 
between AGN and SF luminosity of L1.4 =~ 10 24 5 WHz -1 , and 
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Figure 14. (a) Redshift versus spectral index for those sources of our MOST ATLAS sample, and initial CSS sample with spectroscopic and photometric 
redshifts. (b) Luminosity versus ai?'| for our samples at 1.4 GHz. The pale grey box is the area where SFGs would dominate the population on the basis of 
having a synchrotron spectral index, and a radio luminosity consistent with a star formation rate up to ~ 1000 Mq yr _1 . Any sources in the remaining area 
are AGN. 



Table 6. Properties of our faint initial CSS Sample 





Spectral Index 


Spectroscopic Redshift 


Photometric Redshift 


R magnitude 


1 .4 GHz flux density 


2.3 GHz flux density 


Mean 


-1.58 


0.37 


1.21 


19.03 


1.68 


1.08 


Median 


-1.48 


0.32 


0.97 


18.8 


0.52 


0.36 



! 1.4GHz M08 Lun 



(a] MOST ATLAS Luminosities 



2 



EL 



(b) CSS Sample Luminosities 



Log lt ,(Luminosity (WHz- 1 )) 

Figure 15. (a) Distribution of luminosities for our MOST ATLAS sources 
and (b) initial CSS sample, at 1.4 GHz and 2.3 GHz. 



a spectral index bet ween —0.9 < ot\ ;| < —0.5, typical for SFGs 
dAfonso et al1l2005h . This is a conservative upper limit on the star 
formation rate, assuming that the star formation rate implied from 
the radio emission is very high. The single inverted spectrum source 
with a spectroscopic redshift of z = 1.33 is classified spectroscop- 
ically as a quasar, with broad emission lines. It has an unusual radio 
spectrum that appears to contain an upturn, with Qo;g43 = —1.03, 
and a\X = 1.18, one of only 6 MOST ATLAS sources with an 
upturning radio spectrum. 

The rest-frame 843 MHz, 1.4 GHz, and 2.3 GHz luminosities 
were calculated for each source in the MOST ATLAS catalogue 
with a redshift (spectroscopic or photometric), and the 1.4 and 



2.3 GHz luminosities are shown in Figure [T5h. The rest-frame 1.4 
and 2.3 GHz luminosities were also calculated for the CSS sam- 
ple, and are shown in Figure |15b . The distribution of luminosi- 
ties, in conjunction with their spectral index, is consistent with 
these sources being AGN, as expected. Although it appears in 
Figures [74b and 1 15b that the CSS sample is on average a lower- 
luminosity sample, this is due to the flux limit of the 843 MHz 
MOST data, that preferentially selects the higher-luminosity ob- 
jects. However, as we only find 4 of these MOST ATLAS sources 
in our CSS candidate sample, we are selecting a low-luminosity 
sample of young AGN that is ideal for comparison to our bright 
sample dRandall et alj201lh . 

The MOST ATLAS and initial CSS sample sources with a 
spectroscopic reds hift have also been spectroscopically classified 
as AGN or SFGs jMao et alJliii prep.h . Of the 22 MOST ATLAS 
sources with spectroscopic redshifts, 15 sources were classified as 
AGN, four were classified as SFGs and the remaining two sources 
were unclassified. Three of the SFG sources are discussed further in 
&I7.3I Within the CSS sample, 13 sources were classified as AGN, 5 
as SFGs, and four were unclassified. The location of the CSS sam- 
ple SFGs are shown on Figures [J_6] and [T7] by thick black circles 
surrounding each source. 



7.2 Radio-Far Infrared Correlation 

The radio-far infrared correlation (RFIRC) is a well known tight 
correlation between the far and mid-infrared and radio emission 
from galaxies (primarily SFGs), first noted in the 1970s, and con- 
firme d in later years (e.g. [ Dickey & SaltpeterJI 19841 ; Ide Jong et al.l 
ll985h . lAppleton e t al. (20040 confirmed the correlation holds to a 
redshift 2=1, and it has been recently noted that the RFRIC 
does not evolve above this redshift, but remains relatively con- 
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Figure 16. The radio-far infrared correlation for our sample, and the ATLAS data for CDFS and ELAIS. The N06 and M08 1.4 GHz fluxes are used for all 
sources for consistency. The initial sample of faint CSS sources are also shown, split into candidates with and without redshifts. Objects with a thick black 
circle around the data point are from the CSS sample and are those spectroscopically classified as SFGs. 



stant to a redshift of z ~ 2 dMao et al]|201lh . We have utilized 
the data for both the ELAIS and CDFS field from ATLAS, to 
further investigate the properties of the sources detected in our 
MOST observations, and the initial CSS sample (Figure [T6l>. The 
solid black line on Figure [16] give s the relationship betw een ra- 
dio and infrared flux at 24/^m from I Appleton et alj ( 120041) . where 
924 = log(S24 iim/ S20 cm) ~ 0.84, and the dashed line is the line 
used by M08 to classify a source as unambiguously AGN (above 
the dashed line represented by 924 M m < —0. 16) or a SFG, AGN or 
composite object (below the dashed line). Most of our 843 MHz se- 
lected sources are AGN by the M08 classification. The MOST AT- 
LAS sources which fall below this AGN line are further discussed 
in 3731 

The initial CSS sample is also shown in Figure [TJj] While 
these are clearly low-luminosity AGN on the basis of their spec- 
tral indices (Figure \\Ab), they equally clearly fall predominantly 
in the region occupied by SFGs on the RFIRC. This is an inter- 
esting preliminary result which will be investigated further in a 
follow-up analysis (Randal l et alJlin prepj) . It may be possible that 
the radio spectral index is steep, betraying the presence of an AGN, 
but that the AGN contribution to the 1 .4 GHz flux density is small 
enough (or perhaps is balanced by an equivalent contribution from 
the AGN to the FIR) that the RFIRC is still obeyed. If this is the 
case, the location of the CSS sources close to the RFIRC may 
actually be a consequence of the 1 .4 GHz luminosity being dom- 



inated by star formation, and the host galaxy being a composite 
object. This would be consistent with earlier work suggesting that 
CSS s ources are actively star forming dO'Dealll998l : lLabiano et al.l 
120081 : iMorganti et alj|2009l ; lHoldl2009l) and highlight this popu- 
lation as an ideal resource for exploring AGN feedback effects 
on star formation in galaxies. Only the use of radio discriminants 
(such as morphology or spectral index) will indicate whether the 
objec t hosts an AGN dRov et alJl998l : lNorris. Middelberg & Bovlel 
12007b . 

In Figure [TTJ the 524 values are shown against redshift, 
with s everal diff e rent s p ectral energy distribution (SEP ) track s 
fro m lElviset alJ dl994l) ; iDevriendt. Guiderdoni & Sadatl dl999h . 
and lRieke et alj J2009h overlaid. A sug g ested dividing line between 
AGN and SFGs from ISevmour et alj (2008) is also shown. The 
CSS sample here are again consistent with having their radio and 
FIR emission dominated by star formation as are the three nearby 
MOST ATLAS sources (as labelled in Figures [16] and [TTJ- The 
SED tracks near these objects include those for Luminous Infrared 
Galaxies (LIRGs), Ultra-Luminious Infrared Galaxies (ULIRGs), 
Arp 220 and Mrk 273, two nearby ULI RGs that both host an AGN 
(Engel et af feoi ll ; llwasawa etalj|201ll) . The three MOST ATLAS 
sources in the AGN-dominated region at the bottom of the figure, 
are like ly AGN, and this is supported by the SED of the radio-loud 
quasar dElvis et all 19941) closest to these AGN. Although the SED 
tracks do not extend to the lowest redshifts of our sample, four 
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of the CSS sample are likely to be LIRGs or ULIRGs, similar to 
F00183— 7111, a powerful ULRIG, that has strong ra dio emission 
from both an AGN and vigorous star for ming activity dNorris et al.l 
l201lh . but less luminous in the radio. iNorris et alj j201 lh have 
suggested that F001 83-71 11 is in the earliest ph ase of the forma- 
tion of a quasar, w here a "quasar-mode" AGN (Bes t et al]|2006l : 
ICroton et alj|2006h is hosted by a star forming ULIRG, where the 
star formation is fueled by gas from a past major or minor merger 
event. Whilst the radio jets of F00 183—7111 are currently confined 
by the host galaxy, when they break through this dense gas, the 
jets will begin the transition of this source into a typic al radio-loud 
quasa r, quenching the star formation at the same time 1 Nor ris et al.l 
l201ll) . These CSS sources w ill be further investigated in a future 
paper. iRandall et all din prepj) . 



7.3 Unusual Sources 

We have identified several interesting sources in our final MOST 
ATLAS catalogue, particularly those sources which lie close to the 
radio-far infrared correlation (and consequently are likely to host, 
or even be dominated by, star formation), and one possible ultra- 
steep spectrum object. 



7. 3. 1 Possible SFGs in ATLAS 

Of the four single sources that fall close to the RFIRC, three are 
disk or spiral-like galaxies based on their luminosity profiles in the 
3.6/im SWIRE images (labelled by name in Figures [T6j and!17t. 
which suggests they may be SFGs. The fourth object appears as an 
elliptical, with no evidence of a disk or interaction, and was not fur- 
ther investigated. Figure [T"8l shows the SWIRE 3.6/im images with 
the 1 .4 GHz radio contours overlaid, indicating the clear disk or spi- 
ral structure of these three galaxies. Of the three blended sources 
near the RFIRC, there is less evidence for disks or star-formation 
in the SWIRE data. For the two closest to the RFIRC (S259/S269 
and S1235/S1228/S1230/S1243), the radio emission could possi- 
bly still be associated with star formation, but it is harder to con- 
strain without clear evidence of a disk or spiral structure in the host 
galaxy. We briefly discuss S707, S897 and SI 160 below. 



7.3.2 S707 (ATELAIS J003828.02-433847.2) 

S707 is a possible star-forming galaxy at a redshift of z = 0.048, 
where the host galaxy has a disk structure, as evidenced by is lu- 
minosity profile as seen in Figure [T8h. The radio morphology of 
this object is a single point source at the centre of the SWIRE 
galaxy, suggesting the low-luminosity AGN in the centre of the 
galaxy is the primary source of the radio emission, rather than star 
formation in the disk. Nuclear star formation may also be possi- 
ble, and the optical spectrum of this object suggests there is ongo- 
ing star formation. The radio source has an overall steep spectrum, 
a,fit — —1.17, although it is much flatter between 1.4and2.3GHz 



/„2.3 
("1.4 



-0.45), than 843 MHz and 1.4GHz (aa| 43 = -1.66), 



as expected from an AGN. S707 has a rest-frame luminosity of 
L1.5 = 3.85 x 10 22 WHz -1 , from the observed convolved 1.4GHz 
flux density, resulting in an upper limit on the star formation rate 
of 21M ^yr~ (assum ing all the radio emission is from star for- 
mation; Hopkins 2004). It also has 524 p = 0.57, suggesting a 
small fraction of the radio emission is from star formation, whilst 
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Figure 17. q2ii_im as a function of redshift for the entire ATLAS cata- 
logue, split into spectroscopic and photometric redshifts. The initial CSS 
sample is overlaid as the large red circles for those sources with spectro- 
scopic redshifts, and blue squares for those with photometric redshifts. The 
MOST ATLAS sources with 24/tm detections are shown as the hollow tri- 
angles (sources with spectroscopic redshifts) or inverted hollow triangles 
(sources with photometric redshifts. The SED tracks describe several dif- 
ferent models. The heavy dashed lin e represents an Arp220 type object 
iDevriendt. Guiderdoni & Sadaj|l999h. the thin solid black line is a model 
for Mrk273 iDevriendt, Guiderdoni & Sadatll 19991) . the regular dashed line 
is a radio-loud QS O IeMs et alJI 19941) . the dotted line is a LIRG model 
I Rieke et al. 2009), and the dash-dotted line represents a ULIRG model 
IRieke et al. 2009). The thick black line s hows the dividing line between 
AGN and SFGs from lSevmour et alj J2008h . Objects with a thick black cir- 
cle around the data point are from the CSS sample, and are those spectro- 
scopically classified as SFGs. 



the majority is from the AGN. Optically, this source has a simi- 
lar morphology to the 3.6/im galaxy, with magnitudes B = 15.8, 
V = 15.2 and R = 14.7 from M08. 



7.3.3 S896 (ATELAIS J003429.33-432614.4) 

S896 appears to be hosted by an edge-on disk galaxy, with some 
possible ongoing star formation that can be seen in Figure 1 18b . 
at a redshift of z = 0.053. Although the radio emission ap- 
pears to be emitted primarily from the central object, there is most 
likely some contribution to the overall flux density from star for- 
mation in the disk, as the optical spectrum suggests there is some 
star formation ongoing. This object has a rest-frame luminosity of 
L1.5 = 1.92 x 10 22 WHz _1 , calculated from the observed con- 
volved 1.4 GHz flux density. This implies an upper limit on the star 
formation rate of ~ llMQyr -1 , assuming all the radio emission 
is from star formation. However, it has an overall steep spectral in- 
dex of a fu — —1.21, which suggests that the source of the radio 
emission is a low-luminosity AGN. For this source qi^m = 0.87, 
consistent with the idea that the radio emission is primarily from 
the AGN. 



7.3.4 SI 160 (ATELAIS J003915.1 1-430428.5) 

SI 160 is hosted by a star-forming s piral galaxy, shown i n Fig- 
ureQU:, at a redshift of z = 0.0135 dJones et al.ir2004l2009l) . The 
radio emission is clearly associated with star formation, as the radio 
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Figure 18. (a) 1.4GHz radio contours overlaid on the SWIRE 3.6/im greyscale image for source S707. The radio contours begin at ().4mJy, and increase in 
increments of 0.4 mjy to the maximum level of 2mjy. The cross is the radio position from M08. (b) 1.4 GHz radio contours overlaid on the SWIRE 3.6/jm 
greyscale image for source S896. The radio contours are between 0. 1 and 1 mjy, increasing in increments of 0.2 mjy. The cross is the radio position from M08. 
(c) 1.4 GHz radio contours overlaid on the SWIRE 3.6/xm greyscale image for source SI 160. The radio contour levels begin at 0.5 mjy, and continue up to 
1 mjy, in increments of 0.1 mjy. The large cross is the radio position from M08, and the small cross in the centre of the galaxy is the SWIRE position. 



contours trace the star formation in the spiral arms of this galaxy. 
However, the radio position of the source is offset from the galaxy 
centre, assigned as the location of the brightest peak in the com- 
plex radio structure (M08). The radio luminosity of this source is 
L1.42 = 4.60 x 10 21 WHz -1 which results in an upper limit on 
the star formation rate of 2.5MQyr _1 (if all radio emission is due 
to star formation). The optical morphology of this galaxy is very 
similar to the SWIRE morphology, with a clear spiral structure. 
The optical host galaxy has magnitudes from the SuperCOSMOS 
Science Archive, although these magnitudes should be treated with 
caution, as it is a large, relatively nearby galaxy and measuring ac- 
curate magnitudes in an automated fashion for such sources is dif- 
ficult due to their complex structure. The SuperCOSMOS magni- 
tudes are / = 13.4, R = 12.9 and Bj = 7.9. The optical spectrum 
is also indicative of star formation. SI 160 is detected in all five of 
the SWIRE bands, with a §24^ m = 0.25. This 524 P m value is still 
consistent with the source being either a SFG or an AGN, and it 
most likely there is a non-negligible contribution to the radio emis- 
sion from star formation in this source. The radio flux of this ob- 
ject is interesting, as it has a similar flux density between 843 MHz 
and 1.4GHz (So. 843 = 12.7mJy compared to Si. 4 = 11.6mJy 
or Sa/08 = 12.9 mjy), but drops rapidly away to S2.3 = 5.3 mjy. 
This gives a spectral index ctl'g 4s = —0.2 and a\X — —1-6, with 
an overall spectral index of aju = —1.03. This spectral shape 
is more indicative of AGN activity, but given the radio morphology 
and the FIR constraints, the radio emission from this object is likely 
a combination of star formation and AGN activity. 

7.3.5 S1256 (ATELAIS J003053.26-425215.3): A candidate 
Ultra-Steep Spectrum source 

High redshift radio galaxies are often found by identifying radio 
sources with steep or ultra s teep radio spectra (a < —1), via 
the well-known z-a relation dde Breuck et al]|2004l ; iKlamer et al.l 
1 20061 : Ushwara-Chandra et al j|2010t and references therein;). This 
trend was first noticed in the early 1980s, when it was found that 
the fraction of radio sources optically identified on optical plates 
tended to have the steepest radio spectra. S 1256 is an isolated point 
source, with So. 843 = 4.38 mjy and Si. 4 = 0.81 mJy, that is 
not detected in the 2.3 GHz ATLAS image, with a spectral index 
Q a843 ~ —3.33. Following Zll, using the 3<r detection limit 
of 300/iJy, an upper limit on the spectral index is found to be 
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Figure 19. 1.4 GHz radio contours overlaid on the SWIRE 3.6/jm greyscale 
image for source S1256. The radio contours are at levels of 0.5, 0.75 and 
lmJy. The cross is the radio position from M08. 



ai',4 < —2.41, which suggests that this source is possibly an ultra- 
steep spectrum source, and therefore at high redshift. Further radio 
imaging and flux density measurements would be necessary to con- 
firm S1256 as an ultra-steep spectrum source. S1256 has no optical 
counterpart or confirmed redshift, but is detected in the SWIRE 
3.6 and 4.5^im bands. Its flux densities in these bands are low 
(S3.6 M m = 17.5 /iJy and S4.5/jm = 19.9 ^dy). The infrared mor- 
phology is point-like and is shown in Figure [T9] Another possible 
cause of this steep spectrum is variability; where the radio source 
has faded away since the 1.4 GHz observations were completed, 
and before the 2.3 GHz observations began. Further observations 
would also be necessary to confirm whether variability is the cause 
of the ultra-steep radio spectrum. 



8 CONCLUSIONS 

We have presented a new catalogue of 843 MHz radio sources, 
cross-matched to ATLAS at 1.4 and 2.3 GHz, and explored the 
properties of this catalogue with spectral index as a function of flux 
density. Our results do not support the hypothesis that there is a sig- 
nificant flattening of the spectral index with decreasing flux density 
values. However, we cannot rule out the possibility without further, 
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deep radio data at different radio frequencies. Our analysis of the 
distribution of steep and flat spectrum sources with redshift, lumi- 
nosity and infrared flux density indicates that most of the sources 
in our MOST ATLAS catalogue are AGN, inferred from the spec- 
tral index and infrared properties. An initial sample of faint CSS 
sources in ATLAS has also been selected, with their basic proper- 
ties explored. We will explore the properties of the CSS sample and 
the MOST ATLAS selected GPS sources in depth in future work, 
and compare and co ntrast this faint sample to our bright sample 
(Randall et al. 2011). The spectral index properties of ATLAS ra- 
dio sources across a wide range of frequencies will also be explored 
in future work, with the aim of distinguishing between the propor- 
tions of the populations of low-luminosity, and/or core-dominated 
AGN and SFGs driving this effect. 
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SUPPORTING INFORMATION 

An additional table of Supporting Information is available in the 
online version of this article. 

Table [Q The full version of the catalogue, containing all flux den- 
sity information, ATLAS cross-IDs, optical magnitudes, SWIRE 
IDs and other relevant information. 

Please note: Wiley-Blackwell are not responsible for the content or 
functionality of any supporting material supplied by the authors. 
Any queries (other than missing material) should be directed to the 
corresponding author for the article. 

This paper has been typeset from a Tr-X/ LTr-X file prepared by the 
author. 



© 201 1 RAS, MNRAS OOO.mmi 



